arXiv:1504.04613vl [astro-ph.HE] 17 Apr 2015 


Version 8 : 2015.04.17 


Characterization of the Inner Knot of the Crab: The Site of the 

Gamma-ray Flares? 

Alexander Rudy^, Dieter Horns^, Andrea DeLuca^, Jeffery Kolodziejczak'^, Allyn Tennant‘S, 
Yajie Yuan^, Rolf Buehler^, Jonathon Arons®, Roger Blandford^, Patrizia Caraveo^, 
Enrico Costa®, Stephan Funk^, Elizabeth Hays^^, Andrei Lobanov^, Claire Max^, 
Michael Mayer®, Roberto Mignani^®, Stephen L. O’DelR, Roger Romani^, Marco Tavani®, 

Martin C. WeisskopR 

ABSTRACT 

One of the most intriguing results from the 7 -ray instruments currently in or¬ 
bit has been the detection of powerful flares from the Crab Nebula. Such events, 
with a cadence of about one per year, can be so dramatic as to make the sys¬ 
tem the brightest source in the 7 -ray sky as occurred in 2011 April. These flares 
challenge our understanding of pulsar wind nebulae and models for particle accel¬ 
eration. To pinpoint the production site(s) within the Nebula, a multiwavelength 
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campaign has been carried out using Keck, Hubble Space Telescope (HST), and 
Chandra X-ray Observatory. As the short time scales of the flares (< 1 day) 
suggest a small emitting region, the Crab’s inner knot, located within a fraction 
of an arcsecond from the pulsar, has been proposed to be just such a site. In this 
paper, we focus on IR, optical, and X-ray observations of this feature to see if 
it might be related to the 7 -ray flares and to try to understand its nature. We 
hnd that the knot’s radial size, tangential size, peak flux, and the ratio of the 
knot flux to that of the pulsar as measured with HST are all correlated with the 
projected distance of the knot from the pulsar. A new approach, using singular 
value decomposition for analyzing time series of images, was introduced yielding 
results consistent with the more traditional methods while some uncertainties 
were substantially reduced. 

We exploit the new characterization of the Crab’s inner knot to discuss con¬ 
straints on standard shock-model parameters that may be inferred from our ob¬ 
servations assuming the inner knot lies near to the shocked surface. These include 
inferences as to wind magnetization a, shock shape parameters such as incident 
angle and poloidal radius of curvature Rc, as well as the IR/optical emitting 
particle enthalpy fraction. We hnd that while the standard shock model gives 
good agreement with observation in many respects, there remain two puzzles: 
(a) the observed angular size of the knot relative to the pulsar-knot separation is 
measured to be much smaller than expected; (b) the variable, yet high degree of 
polarization reported is difficult to reconcile with a highly relativistic downstream 
how. In contrast, the IR/optical hux of the inner knot is marginally consistent 
with the scenario that the shock accelerates most of the optical emitting particles 
in the nebula. 

Subject headings: neutron stars: general — pulsars: individual — Crab nebula 
and pulsar (PSR B0531-I-21) — inner knot 


Introduction 


The Crab pulsar and its nebula are one of the most studied targets in the sky at 
all wavelengths. It serves as a test bed for pulsar theories as well as, more generally, for 
astrophysical non-thermal processes. A very thorough review of the Crab system, including 
high spatial resolution observations with HST and Chandra, showing the complex, d ynamic 
i ntera ction of the pulsar wind wit h the surrounding rn e dium , has been compiled by [Hester 
(120081) . A more recent review by iBiihler fc BlandfordI (120141) also includes the saga of the 
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7 -ray flares. 


A crucial question implied by the observations of the 7 -ray flares is “how are the parti¬ 
cles being accelerated to sustain the flaring behavior?” Radiating particles must have PeV 
energies to emit synchrotron radiation at the observed 7 -ray wavelengths. A mechanism is 
needed that can accelerate particles to such high energies in less than one gyration time in 
the ambient magnetic held. For example, the large hare of 2011 April produced a radiant 
energy equivalent to the energy stored in the mean magnetic held within a region of order 
2 X 10^'^ m subtending an angle ~ 0.3" at 2 kpc, consistent with the emission-region size 
estimated from light-travel-time arguments fjWeisskopf et al.ll2013l) . 


An intens e theoretical ehqrt is currently being devoted to explain the variab 
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20131 1. However, the mechanism driving the hares, their impulsive nature, the ~ 12-month 
recurrence time, and the location, remain unknown. Possible explanations include plasma 
instabilities in the nebula, magnetic reconnection, discontinuity in the pulsar wind acceler¬ 
ation, emission from the anvil, emissio n from the inn er knot, or portions of the termination 
shock. See, for example, the review by iHesterl (120081) for the association of specihc features 
with this nomenclature 

There are a number of reasons to suspect that the feature kn own as the ’’inner knot 


shoul d be considered as a possible site for the origin of the hares fe.g.. lKomissarov fc Lvutikov 
2OIII ). First, the luminosity of the hares accounts for a relatively high fraction of the pulsar 
spin-down luminosity (^ass uming a radiation efficiency of 3-5%), suggesting an origin in a 
region close to the pulsar (ITavani et al.l 1201 111 . Second, the very short variability time scales 
combined with causality arguments, stro ngly constrain the size of the hare’s emitting region 
( Buehler et al. 2012 : Striani et al. 2011 1. 


A multiwavelength campaign to study the Crab, described in IWeisskopf et al.l (120131) . 
is being performed. Another 7 -r ay hare—the secon d largest observed to date—was detected 


with Fermi/LAT in 2013 March (iMaver et al.ll2013l ). This triggered further observations with 
Keck, with the Hubble Space Telescope (HST), and with the Chandra X-ray Observatory. 
Here we concentrate on what was learned about the inner knot as a consequence of this 
campai gn. We note that very-high-energy (E > 100 GeV) co ntemporaneous d ata from the 
HESS ( H. E. S. S. Collaboration et al. 2014) and VERITAS ( Alin et al. 2014) Cherenkov- 
telescope arrays found no signihcant correlation of TeV hux with the hare observed with 
Fermi/LAT. 
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The observations are presented in ^ The measured properties of the knot are given 
in ^ §1] compares observations of the inner knot at the different wavelengths. Comparison 
and examination of possible correlations of the knot’s properties with 7 -rays detected using 
the Fermi/LAT are discussed in §3 §S] discusses theoretical implications, ^summarizes. 


The Observations 


The feature we refer to a s the “inner knot” (aka “knot 1”, or “synchrotron knot”) 


was discovered by iHester et all fjl995[) in a detailed study of the Crab using high-resolution 
images taken with HST. The feature, although resolved in those data, is small (~ 0.1") 
and close to the pulsar 0.6"). The feature was hrst suggested to be associated with 
a shock or instability in the jet ~ 1500 AU from the p ulsar, due to its good alignment 
with the jet and elongation perpendicular to the jet. Later iKomissarov fc Lvubarskvl (120031) 
proposed, based on two-dimensional magnetohydrodynamic (MHD) simulation of the pulsar 
wind nebula, that the knot represents Doppler-boosted emission fro m an oblique termination 
shock. Synthetic synchrotron images fr om these 2 D simu lations (IKomissarov fc Lvutikov 


20 111) and more recently 3D simulations (jPorth et ahl 120141) indeed show a knot-like feature 
when the system is viewed from the proper (and reasonable) orientation. Moreover, the 
observed high degr ee of polarization and a position angle aligned with the symmetry axis 
(IMoran et al.l 120131) seems to lend further support to the oblique-shock scenario. 


The salient features of the inner knot deduced prior to our observati ons are tha t it ha s 
a power-law spectru m oc z/“" with index a ~ 0.8 in t he IR to optical (ISollermanI (120031) . 
Melatos et al.l (l2005l) j and, as most recently measured bv iMoran et al.l (120131) . the flux varies, 
the position appears to change, the degree of polarization is very large but does not appear 
to change with the flux (but see the discussion at the end of j|5]), and the magnetic held 
(determined from the position angle of the optical polarization) is orthogonal to the axis of 
symmetry (presumably the spin axis) and consistent with the direction of the held for the 
rest of the inner nebula. 


2.1. Keck 


We obtained K’ and H band near-infrared (NIR) images of the inner 40" of the Crab 
Nebula using the Keck Ne ar-Infrared Camera 2 (NI RC2; K. Matthews, PI), with laser-guide- 


star adaptive optics (AO: IWizinowich et ah II2006I) . From 2012 February to 2014 January, 
we obtained 11 observations of the nebula during 7 -ray quiescent periods and 1 target- 
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of-opportunity observation (K 8 ) triggered by 7 -ray flares. Table [T] summarizes the Keck 
observations. 


Data reduction was carried out using the pipeline developed bv iGhez et al.l (120081) . In¬ 
dividual frames were sky-subtracted, flat-fielded, bad-pixel-corre cted, and co r rected for the 
NIRC2’s camera distortion. We used the solution developed bv lYelda et al.l (I 2 OIOI) . which 
matched the positions of stars in the globular cluster M92 measured with HST to those 
measured in the NIRC2. Frames were checked to ensure that they showed no irregularities 
in the point-spread function (PSF), and that the Strehl ratio was at least ^ O-IO- In di- 
vidual frames were then combined using the Drizzle algorithm flFruchter fc Hookll2002l) to 
provide sub-pixel accuracy and high signal-to-noise images for analysis and modeling (1 13.11 
and Appendix). 


The observations were conducted with two different pixel scales, using the NIRC2 “wide” 
camera, with 0.04" pixels and a 40" x 40" held of view, and the NIRC2 “narrow” camera, 
with 0.01" pixels and a 10" x 10" held of view. The total effective exposure times ranged 
from 420 to ~ 2600 seconds. The shorter exposures result in optimal use of 1-hour Target- 
of-Opportunity interrupts—including acquisition, calibration, and sky measurements. 


The spatial resolution of Keck’s NIRC2 with AO proved to be very effective in resolving 
the detailed structure and position of the knot. Figured] shows the 12 NIRC2 images. In the 
narrow camera, the center of the pulsar and the center of the knot are separated by ~ 60 
pixels. Typical values for the signal-to-noise ratio of the knot are 5. The PSF of NIRC2 with 
AO varies slightly in time due to the variable performance of the AO system. As the pulsar 
and knot are well resolved from each other in the Keck images, we chose not to deconvolve 
the PSF from the Keck images. 


2.2. HST 

Since 2012 January, we have obtained 17 observations using the HST Advanced Camera 
for Surveys (ACS) Wide-Field Channel (WFC), with 0.05" pixels. Each observation com¬ 
prises 4 exposures totaling 2000 seconds and uses a standard 4-point BOX dithering pattern 
to All the inter-chip gap as well as to allow efficient cosmic-ray cleaning. We employed the 
F550M filter, which is well suited to sampling continuum emission from the Crab with almost 
no contamination from line emission. 

Table [2] summarizes the observations, listing the date and the mean observing time. The 
Appendix discusses in detail our analyses of the 17 HST/ACS images. Figured displays the 
processed images, which have removed effects of the HST/ACS WFC point spread function 
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using singular-value decomposition (SVD) and Richardson-Lucy deconvolution. 


2.3. Chandra 


Our approach for searching for X-ray emission from the inner knot with Chandra in¬ 
volved a special use of both the high time resolution and high spatial resolution of the High- 
Resolution Camera for Spectroscopy (HRC-S) array in timing mode. So doing combines the 
best spatial resolution and the best time resolution afforded by the Chandra instrumenta¬ 
tion. To achieve the time resolution, it is vital that the counting rate be below the rate 
where telemetry is saturated (see, e.g.. iTennant et al.ll200ll) . Consequently, the Low-Energy 
Transmission Grating (LETG) was inserted into the optical path to reduce the flux. As in¬ 
serting the grating is of itself insufficient to avoid telemetry saturation, the trigger threshold 
of the HRG-S was was raised to reduce the rate of detected events. This approach was hrst 
successfully used in a trial experiment (ObsID 11245) on 2010 November 16. 

The next set of observations in this mode was initiate d by us on 2013 March 5 as part 
of a sequence triggered by a y-ray flare fjOiha et al.ll2013l) . Unfortunately, the gain of the 
HRG had dropped since 2010. Gonsequently, the trigger threshold was too high during 
our hrst observation on March 5-6, reducing the counting rate from the pulsar to about a 
quarter of that expected, dramatically decreasing the number of detected counts. The next 
HRG observation was to take place on March 10 and the observation had already been fully 
planned with the commands to execute onboard the spacecraft. To replan at this stage 
would require an extremely short turnaround, yet it was accomplished. Finally, beginning 
on 2013 October 22, we triggered ano ther sequence of H RG observations in response to the 
announcement of another y-ray hare (IBuson et al.l 120131) . 


We also consider ObsID 9765 taken on 2008 January 22. For this observation the HRG 
was operated in a more standard conhguration with the threshold set at the nominal value. 
However, to reduce the count rate, a blade was inserted into the optical path, which is no 
longer allowed. This blade also reduced the dihracted hux in one direction. By chance, this 
reduction is greatest near the nominal location of the knot. This reduction, together with a 
long integration time for this observation, results in ObsID 9765 providing the best data for 
study of the knot. Table [3] summarizes the relevant Chandra observations. 
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3. Properties of the Knot 
3.1. Keck 

Table 0] lists the best-fit values of the model parameters for the inner knot for each of 
the 12 Keck observations based upon the analysis and modeling described in the Appendix. 
The properties are these: 'ipOi the position angle (E of N) from the pulsar to the center of 
the knot; ro, the projected separation between the pulsar and the knot; FWHMr, the full 
width at half maximum of the knot in the radial direction; FWHMt, the FWHM of the knot 
in the tangential direction; and Fk/Fp, the ratio of the flux in the knot to that in the pulsar. 
The flux ratio is used as measuring the fluxes of the 5 stars in in the NIRC2-wide camera’s 
held of view shows that those huxes change relative to each other from one observation to 
the next. This occurs because the nebular background varies spatially across the detector 
in a non-trivial way, making it difhcult to determine precisely the hux of any given source. 
Consequently, we measure the ratio of the knot’s hux to the pulsar’s, rather than the knot’s 
absolute hux. Figure |3] plots these properties versus the pulsar-knot separation, based upon 
the Keck data. 

We seek possible correlations of these properties with the separation tq from the pulsar, 
of the following form: p(ro) = p(ri)[ro/ri]'^. We choose a reference distance ri = 0.6692", 
obtained from the logarithmic average of the 12 values of tq. The position angle of the knot, 
'ipo, was ht to the linear regression = f/^o(^i) + V’o(^i)[^o ~ ''"i]- Table [5] tabulates the 
results of the regression analyses of the inner knots intrinsic properties with its (projected) 
separation from the pulsar. Based upon the F-test probability, the power-law regression 
analysis hnds marginally signihcant correlations of the radial width, tangential width, and 
peak intensity (surface brightness) with separation. The correlation of the knot’s flux to the 
separation is not signihcant. 


3.2. HST 

Table [6] lists the best-ht values of the model parameters for the inner knot for each of 
the 17 HST observations based upon the modeling and analysis described in the Appendix. 
The properties are these: the position angle i/jq of the knot with respect to the pulsar, 
the (projected) distance tq from the pulsar, radial width (FWHM^) and tangential width 
(FWHMt) of the knot, the peak surface brightness Sk, the knot hux F^, the pulsar hux, Fp, 
and hux ratio F^/Fp. Although we expected the pulsar hux to be steady, we measured a 
small {2.5a) variation. Consequently, we prefer to work with the hux ratio Fk/Fp. 


In the Appendix we discuss two ways of analyzing these data. Table [6] presents the 
results for one of these methods. However, our general conclusions concerning the knot’s 
optical properties are for the most part the same for either analysis. Figure H] plots the knot 
properties versus the pulsar-knot separation, based upon the HST data. 


As with the Keck data, we perform a regression analysis of the inner knots intrinsic 
properties with its (projected) separation from the pulsar (Table [7]). We choose a reference 
distance ri = 0.6559", obtained from the logarithmic average of the 17 measured values of vq. 
The position angle of the knot, was £t to a linear regression 'ipQ = +t/^o(^i)[^o 


Based upon the F-test probability, the power-law regression analysis finds statistically 
signihcant correlations of the t angential and radia l widths with separation from the pulsar, as 
well as the previously known flMoran et al.l 120131) anti-correlations of flux and peak surface 
brightness with separation. The correlation of i/jq with separation vq is the least strong. 
There is a small systematic bias in the power-law index of the FWHM^ correlation, in that 
the radial width of the knot at times approaches the resolution limit of the instrument and 
analysis procedures. Applying our analysis procedures (see Appendix) to simulated knots 
of known true radial width, we estimate that the smallest measured values of FWHM^ are 
oversized by about 7% and the largest by about 3%, which would increase the power-law 
index for the FWHMr-ro correlation from q = 0.72 to about 0.8. Finally, we emphasize that 
we hnd similarly strong correlations—except for the flux—using the more traditional HST 
analysis methods described in the Appendix. 


3.3. Chandra 

For the Chandra data, the main question is whether X-ray emission from the knot can 
be detected—especially in the vicinity of the bright. X-ray-emitting pulsar. The answer, 
unfortunately, is “No” and here we set an informative upper limit to the X-ray flux of the 
inner knot relative to that of the pulsar. 

Table E] shows the number of counts in each of three extraction regions: (1) centered 
on the approximate location of the optical knot, (2) centered on the pulsar, and (3) placed 
on a reference region to estimate the background. Details may be found in the Appendix 
and are illustrated in Figured The analysis was performed both for the full phase-averaged 
data and for the data at pulse minimum. Consider, for example, the phase-averaged data 
for ObsID 11245. Comparing the inner knot counts to the background (Table [8]) we hnd 
an excess of 167 counts. If the excess were really due to the knot, we would accordingly 
expect an excess of 50 (= 167 x 0.3) counts at phase minimum. As this was not the case, we 
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attribute the excess to the wings of the PSF from the pulsar and posit that it is due to the 
spring/fall asymmetry discussed in the Appendix. As the impact of that asymmetry is below 
the statistical noise for the data taken at phase minimum, we then use the phase-minimum 
data to derive upper limits. We convert our phase-minimum upper limit into a phase-average 
upper limit by accounting for the 0.3 phase duration of the minimum and then divide by the 
phase average pulsar flux. Results are listed in the last column in Table El 


4. Comparison amongst observations of the inner knot 

Our Chandra observations were unable to detect X-ray flux from the knot, setting a 
3-sigma upper limit to the ratio of knot flux to pulsar flux Fk/Fp < 0.0022. This upper limit 
is 20-40 times smaller than the ratio measured in the optical {Fk/Fp = 0.044 ± 0.001) and 
in the near-IR {Fk/Fp = 0.079 ± 0.005), indicating that the optical-X-ray spectral index of 
the knot is signihcantly steeper (by at least 0.5) than that of the pulsar. 

Both our infrared and optical observations detected time variability in the pulsar-knot 
separation, with roughly consistent separations when measured contemporaneously (see the 
lower portion of Figure [6]). Correlations of knot properties with pulsar-knot separation are 
very strong for the optical measurements but rather weak for the infrared measurements. 
We attribute this difference to a combination of factors—including poorer statistics for the 
much shorter NIR measurements and a number of systematic effects, resulting in part from 
the NIR PSF and background varying across the held and in time. As the lower portion 
of Figure Eland Table [4] indicate, NIR-measured separations at essentially the same epoch 
exhibit a scatter substantially larger than the quoted statistical errors in the measurement, 
thus suggesting the presence of unaccounted systematic errors. 


5. Comparison with y-ray fluxes 


It is clear from the previous sections that the intrinsic properties of the optical knot are 
correlated with its projected separation from the pulsar. The correlation for the NIR knot 
is but marginally signiflcant. Thus, in searching for any relationship between the knot and 
the y-ray flux, we concentrate on the separation versus y- ray flux. The analy sis performed 
to reduce the Fermi/LAT data is the same as described in iMaver et al.l (120131) . 


Table El lists the dates and separations measured in the Keck and in the HST observa¬ 
tions, followed by the corresponding 12-hour-average Fermi/LAT > 100-MeV fluxes. The 
upper portion of Figure Elshows the Fermi/LAT flux as a function of time covering the inter- 
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val when our Keck and HST observations were made. The vertical lines identify times when 
observations were made with Keck (red) or with HST (green). The lower portion of Figure [6] 
plots the pulsar-knot separations on the same time scale. Figure [7] expands Figure [H] at the 
times of the largest y-ray flares. 

Figure [8] plots the y-ray flux versus the measured radial separation of the knot from 
the pulsar. The results captured in this hgure clearly do not represent a signihcant simple 
correlation of the knot properties with the y-ray flux: The formal correlation coefficient is 
only 0.2. 

Inspection of the y-ray light curve and the time series of pulsar-knot separations (Fig¬ 
ure [6]) indicates that the most energetic y-ray flares occurred at a time which is not precisely 
coincident with the time of greatest pulsar-knot separation. It is reasonable to ask how 
probable such a coincidence would be as a chance occurrence. A simple answer is that the 
probability of a chance coincidence is just the difference in the central times of the two events 
divided by the total duration of this observing program, as described below. 

We next estimate the time and duration of a y-ray flare and the excursion in tq (lower 
portion of Figure |6]), using a simple Gaussian model in both cases and including a linear 
baseline trend for the y-ray flare. This is not to suggest that a Gaussian properly describes 
the underlying physical processes, but simply to enable comparison of the two time series, one 
for the flares, one for the separations. Fitting separation data in the MJD range 56180-56400 
and y-ray fluxes in the range 56345-56370 places the midpoint of the inner knot excursion 
at 56331.8, followed by the y-ray flare at 56357.3, 25.5 days later. The Gaussian width is a 
measure of the event duration with FWHMs of 70.7 and 7.6 days for the separation excursion 
and y—ray flare data respectively. 

The time between our hrst and last pulsar-knot separation measurements is 826.1 days. 
However, on 4 occasions the gap between separation measurements was longer than the 
70.7 day duration of the observed separation excursion. We assume then that a separation 
excursion would be observable if an observation fell within the FWHM window and not 
observable if it fell outside the window. To obtain the total observing program under this 
assumption, we remove the time of 4 gaps (140.1 days) that exceeded 70.7 days making the 
total duration of the observing program 686 days. Thus an estimate of the chance probability 
is 25.5/686 = 3.7%, which is approximately at the 2a level. 

It is worth noting that the largest pulsar-knot separation incursion (toward the pulsar) 
occurs near the second most energetic y—ray flare. However, the knot separation data are 
sparser at this time. Therefore, we use the previous FWHM value of 70.7 days to establish 
the time of the incursion. With this prior, htting the separation data in the MJD range 


11 


56390-56670 and 7 -ray fluxes in the range 56575-56590 places the midpoint of the inner knot 
excursion at 56566.0, followed by the 7 -ray flare at 56582.6, 16.6 days later. 

In this second case, the comparable chance probability estimate is 16.6/686 = 2.4% 
(~ 2.5a). As each variation of the knot-pulsar separation is of the opposite sign, we must 
multiply the probabilities by a factor of 2 when combining them. So our estimate of the 
chance probability that the 2 largest knot excursions/incursions would occur within the 
observed intervals of the 2 largest 7 -ray flares is 0.037 x 0.024 x 4 = 0.4%, which is approxi¬ 
mately at the 2.9a level. There are additional 7 -ray flares; however, the third most energetic 
flare occurs only 12 days later than the one on MJD 56582.6 and therefore is well within 
the assumed 70-day range of the associated knot incursion. The other flares are significantly 
weaker and less well-sampled by the optical/infrared observations. 

Although coverage is sparse, there is no compelling evidence in the pulsar-knot sepa¬ 
rations for large (> 0 . 01 ") excursions from the mean separation at time-scales comparable 
to the flare durations. Indeed, knot incursion/excursion time scales appear to be > 10 x 
longer than flare durations, so direct correlation of the 7 -ray light curve with the separation 
time-evolution produces correlation coefficients that are rather difficult to interpret without 
associated simulations. Such simulations would require assumptions regarding an underlying 
physical model relating the two curves. As we have no such model, any correlation analysis 
would be purely ad hoc. 

The 0.4% value for the chance probability is based upon a rather simplistic analysis 
involving only 2 flares and was, of course, determined a posteriori. However, in continuing 
to monitor the Crab at multiple wavelengths with the goal of understanding the origin of 
its 7 -ray flares, the inner knot will continue to be a feature of interest. Based upon the ro 
time scale of 71 days (30-day Gaussian la), a monitoring program which includes monthly 
optical observations is necessary to detect significant inner-knot excursions from its mean 
pulsar separation. Monitoring every ~ 3 weeks should be sufficient to characterize fully these 
variations and to verify the small statistical sample from which these parameter estimates 
were derived. 


We also note that an inner-knot incursion with an 2-month FWHM t imescale and 
0.075" ± 0.025" pulsar-knot separation inward amplitude was also reported bv iMoran et ah 


(l2013f) . This event appears to be quite similar in duration and amplitude to our observations. 
The reference also reports an cs 2a apparent change in the knot polarization parameters, 
which occurred 22 ± 5 days after the time of maximum incursion. The degree of linear 
polarization steps from 59.9% ± 6.9% on 2005 Nov 25 to 42.8% ± 6.2% 2005 Dec 5 and the 
position angle changes from 125.4° ± 3.3° to 121.6° ± 4.2°. Interestingly, the time frame of 
this polarization shift occurred in a window closely matching both the 25.5 and 16.6 day 7 - 
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ray flare delays after the peak inner-knot excursion/incursion times discussed above. While 
the claimed 2a result is not very signihcant, we observe that prior to the step near the start 
of 2005 December, all 10 of the reported measurements of the degree of linear polarization 
fall into a population with mean 60.7% and small standard deviation 1.0%, whereas the 
2 subsequent measurements had a mean of 43.4% with standard deviation of 0.8%. If we 
suppose that the quoted errors in the reference apply to the absolute, but not relative 
measured values, then the size of the step change would be > lOa instead of ~ 2cr. This 
may be a reasonable assumption, as the point-to-point measurement variations are smaller 
by a factor > 6 than the quoted errors. Similar arguments would apply to the polarization 
position angle. From this we conclude that polarization should also be of interest in future 
attempts to relate inner-knot behavior to y-ray flaring. 


6. Implications for Theoretical Models 


In the basic model of the Crab Nebula fe.g.. lRees fc Gunnlll974l) . the pulsar wind passes 
through a shock at a radius ~ 3 x 10^^ cm where the wind momentum flux balances the 


nebular pressure. However, the w ind shoulc 


be a nisotropic. For example, the momentum 


flux in a split-monopole oc sin^ 6 ^ (Bogovalovlll999) and recent simulation s suggest oc sin'^ 6 ' 


flTchekhovskov et al.ll2013l: iTchekhovskov, Philippov, fc Spitkovskv II2015I) where 6 is the po¬ 
lar angle. As a result, the shock is likely to be quite oblate. The sections of the shock near 
the poles would be oblique and much closer to the pulsar than the equatorial part. The 
observed radiation presumably comes from relativistic particles (electrons or positrons) ac¬ 
celerated behind the shock. Since the outflow from an oblique shock can remain relativistic, 
we would be able to see a compact, emitting feature (the inner knot) if the relativistic outflow 
happens to be aligned with our line of sight so that its emission is beamed. The inner knot 
should appear to us as having some offset from the pulsar due to the deflection of the outflow 
from its presumed radial motion (Figure [9]). If the shock is approximately axisymmetric as 
one would expect, the projected emitting site would fall on the projected symmetry axis, 
leading to an alignment with the jet. The shape of the inner knot should also be more or 
less symmetric about the axis; whether it is elongated parallel or perpendicular to the axis 
depends on the geom etry of the oblique shock. Furt hermore the synchrotron emission should 
be linearly polarized (iKomissarov fc Lvutikovll201ll. Yuan & Blandford 2015 in preparation). 


The shock model predicts that the scaling between the observed properties of the knot— 
projected pulsar-knot separation tq, tangential angular width FWHMi, radial angular width 
FWHMr and surface brightness Sk —should be determined by the upstream magnetization 
a = Bl/noni'y^mc^ (where m is electron mass and Bi, ni, 71 are upstream magnetic held. 




















13 


proper density, bulk Lorentz factor, respectively) plus three more parameters that charac¬ 
terize the shape of the shock near the emitting site: The unprojected pulsar-knot separation 
the incident angle (5i, dehned to be the angle between upstream velocity and the shock 
surface as shown in Figure 0 and the shock poloidal radius of curvature -Re- 


Let the outflow be deflected from the radial direction by an angle A at the shock. We 
then have tq = r^A. The knot tangential size is determined by the Doppler beaming and 
the shock radius of curvature in the toroidal direction as FWHM^ 2rk/'-y, where 7 is the 
downstream Lorentz factor. In the simplest shock model with an isotropic plasma, both 
A and 7 are functions of a and (5i: A = (5i — arctan(xtan(5i) and 7 = 1 /(sin (5i-\/l — y^), 
where y = V 2 ±/vi± = Bi/B 2 = (1-|- 2 (T-|- a/IGct^ -|- 16cr -|- l)/6(l-|-cr) is the compression ratio 
at the shock, ui 9 I being the component o f upstream/downstream velocity perpendicular 
to the shock (IKomissarov fc LyntikovI 120111 ). The radial knot size FWHMr should also be 
proportional to but will have an additional factor involving the ratio Rc/^k- 


As to the intensity, we assume that the emitting particles have a power-law distribution 
oc in the fluid rest frame, where n', 7 ' are particle density and Lorentz fa ctor in the 

same reference frame, and p ~ 2.6 from the measured spectral index of the knot fjSollerman 


2 OO 3 I: iMelatos et al.l 1200511 . 


The emissivity in the fluid rest frame is j//(zz') oc n', where B' and 
u' are magnetic held and radiation frequency measured in this reference frame. After 
transforming to the nebula frame, the emissivity has a dependence on the how Lorentz 
factor 7 , particle (improper) number density n, and magnetic held B as follows: oc 

D^+(P“i)/^(n/ 7 )(-B/ 7 )(^’i'i)/^, where V = vjv' = -^/l — (3"^ /(1 — /3 • n) is the Doppler factor, 
in which f3 is the how velocity and n is a unit vector along the direction of line of sight. 
The surface brightness corresponds to the emissivity integrated along the line of sight in 
the nebular frame: = f j^dl. To estimate the peak surface brightness, one takes the 

maximum Doppler factor V 27 , and we get S^^p^ak cc where i is an estimation 

of the length of the emitting region along the line of sight. As roughly B oc and n oc r/"^, 
we have Sy^peak oc 


These previous relations allow us to set the constraints discussed below using both the 
steady and variable properties of the knot. 


6.1. Variability 

In the shock model, variability can arise either upstream or downstream of the shock. 
Firstly, stress tensor variation in the nebula can cause the shock radius and shape to vary 
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over time. Three dimensional MHD simnlations fjPorth et al.l l2014l) show that even if the 
npstream condition is hxed, the post-shock flow can be qnite variable: there is vortex- 
shedding from the termination shock, and the shock constantly interacts with waves and 
vortices in the nebnla. Those variations are expected to be mostly magnetic in origin, and 
the interactions between the shock and t he downstream flow are qnite nonlinear. The pictnre 
is similar to earlier 2D simnlations fe.g.. lCamns et al.ll2009h except that the kink instability 
sets in so that the hoop compression is less prominent as in the 2D case; also short-term 
variability is less prononnced. The 3D simnlations give a typical variation time scale of less 
than one year. 


Now snppose that the change of downstream pressnre canses the shock radins to 
change. For simplicity we hrst assnme that A, i, 7 and p stay more or less the same, then 
from the above scaling relations, we have S^^peak cc FWHMf oc tq, and knot flnx 

oc S'y^peafcFWHMt FWHMr oc assnming FWHMr to be also ronghly proportional 

to rfc. In reality A, i, 7 change with the shock radins as well, as the shock shape and 
npstream incident angle vary. Conseqnently, the observed correlation is more complicated 
and less clean; however, the general trend is consistent with onr expectation of the shock 
model. 


Variations initiated by the pnlsar can also canse the knot properties to change. In snch a 
scenario, the shortest possible variation time scale is F ~ rfc(l — cos(l/ 7 ))/c ~ rk/{2'y‘^c) ~ 
1 day(rfc/10^^cm)( 7 / 5 )“^ np to 10% inaccnracy assnming 7 ^ 3 . However, there is no 
evidence in the pnlsar timing that this is actnally happening. 


6.2. Time averaged properties of the knot 

The properties of the knot may be nsed to set constraints on the flow composition 
and how mnch dissipation occnrs at the shock. One constraint comes from geometrical 
relations amongst the three measnred angles vq, FWHM^, and FWHM^. For the sim¬ 
plest shock model, we have tq = r^A and FWHM^ 'r'k/j, where A and 7 only de¬ 
pend on a and hi, snch that the ratio FWHMt/ro sets limit on a and hi. We hnd that 
FWHMi/ro = 2 /( 7 A) = 2 sin hi V^W-^/[hi — arctan(xtanhi)] > 2 . 8 , with the minimnm at 
a = 0 and hi = 0.39. In contrast, the observation gives FWHMi ~ O.Sro, indicating that 
some important ingredient is missing in the standard shock model. For example, the np¬ 
stream flow might have already been bent before arriving at the shock; and/or the anisotropic 
particle distribntion downstream might be important in determining the size of the knot. 
Althongh more detailed investigations will be described elsewhere, the above discnssion on 
variability and the following discnssion on energetics shonld be valid for qnite general shock 
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models. 


The shape of the knot also has some interesting implications. In both the Hubble and 
Keck images, the knot is not just a simple ellipse that only involves quadratic terms of 
coordinates; it has curvature and appears to be convex away from the pulsar (“smile”), so at 
least a third moment of the coordinates is needed to describe the shape. A simple intuitive 
picture is based on the geometry that the shock surface is similar to a torus centered on 
the pulsar and we look through the hole from the bottom. If emission from the shock 
outflow roughly follows the local toroidal magnetic held, one would imagine the knot to be 
convex toward the pulsar (“frown”)—contradictory to what we see. In fact, we should not 
overlook the possibility that the outhow may have signihcant bending downstream. If the 
outhow bends toward the equator by an angle ~ I /7 within the emission length, we would 
see the knot more extended (contracted) on the far (near) side with respect to the pulsar, 
as additional how lines enter (originally aligned how lines leave) the emission cone on the 
corresponding side. 


Another constraint comes from the observed hux of the knot. In the K' band, the typi- 
cal knot-pulsar hux rat i o is 0.063 (Table 0]). If we take the photometry of the pulsar from 


Sandberg fc SollermanI (120091) . the pulsar Kg ban d magnitude is 13.80 ± O.OI, after dered¬ 


dening with the best ht va lue E{B — V) = 0.52 fjSollerman et al.ll2000l) and the extinction 
law from IFitzpatricki ( 19991) . this corresponds to a hux of 2.2 x I 0 “^®ergs“^cm“^Hz“^. So the 
inner knot hux in K' band should be roughly K^not _ ^ I0“^^ergs“^cm“^Hz“^. Similarly, 


-1 


optical data (Table[H]) on average give a dereddened hux of = I x I 0 “^^ergs“^cm“^Hz 
However, we do not attempt here to get a spectral index using the time averaged huxes due 
to the variability of the knot and the fact that most of th e mult i-wa velength observations 
are not simultaneous. Previous measurements bv ISollermanI (120031) and iMelatos et al.l (120051) 
indicate that the spectrum of the k not in IR to optical range is a power law F^, oc with 
0 . 8 , and IMelatos et al.l (120051) further measured that this power law extends to NUV 


a 


(160 — 320 nm), with a hux in NUV = 0.12 mJy = 1.2 x I0“^'ergs“^cm“^Hz“ —this 
is the highest energy at which the knot has been measured. A few more recent measu re¬ 
men ts show slightly diherent spe ctral index ranging from 0.63 (iTziamtzis et al.l 120091) to 
1.3 (ISandberg fc SollermanI 120091) but unfortunately these are based on non-simultaneous 
data. It would be very interesting to rehne the spectral measurement using near simulta¬ 
neous IR/optical images in the future. In what follows we adopt a hducial spectral model 
F^not ^ I 0 “^^(i//I 0 ^^Hz)“°'®ergs“^cm“^Hz“^ with the caveat that this could change with 


better observations. 


Using the above spectral model, we hnd that the knot integrated hux in IR/optical band 
((O.I — 2) X 10^^ Hz) is Fq/ir ~ ^ 10“^^ ergs“^cm“^. Thus, the luminosity of the knot 
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per steradian is Lo/m,n = -D^Fq°j^ = 10 ^^ergs where H ~ 2 kpc is the distance of the 
Crab. If we dehne e = Lo/iR,r 2 /(i?/ 47 r), where E = 5 x 10^® ergs“^ is the pulsar spin-down 
power, e can be used as a measure of the radiative efficiency of the shock along the direction 
of line of sight in IR to optical range. Here we hnd e = 2.5 x 10“® and the shock is highly 
adiabatic. 


We can get an estimation of the enthalpy fraction of IR/optical emitting particles in 
the downstream flow 

/ n _ 8 7fcooi ^ 8 

3 tfiow E/{47r) 3 fflow E/(47r)’ ^ 

where tflow is the flow time scale and = Gneom^c^/^ years 

is the synchrotron cooling time of the particles in the fluid rest frame (we have adopted 
B = R_ 3 mG and u = for the numerical value). Thus the particle injection rate 

per unit steradian in corresponding energy band is 


7V(7') 


3h(y) E 

A'y'mc^ 47r7 


( 2 ) 


Adopting typical values 7 = 675, fflow = 'ri 7 l 0 ^^ cm/c, we hnd that for IR emitting particles, 
7 ir = 0.057fiir3^^ry, and Nir = 2 x lO^^'y^Bzlr^r s“^sr“h The spectral index of the 
knot indicates that the particle distribution powerlaw index should be p ~ 2 . 6 , thus the 
IR emitting particles should comprise the majority of particle pressure downstream. From 
the above estimations, it seems that the particle injection rate at the shock is marginally 
consistent with the scenario that most of the IR/optical emitting particles are provided by 
the shock. 


Regarding particle acceleration mechanisms, we notice that IR emitting particles usually 
go through N = Larmor orbits, where Ug = vb/z' is the relativistic gyrofrequency— 

the cyclotron frequency {vb) divided by the electron’s Lorentz factor ( 7 '). Eliminating 7 ' 
in favor of the magnetic held {B) and observing frequency (z/), the number of orbits within 
the how timescale is 6 x . To allow sufhcient time for acceleration 

mechanisms requiring instabilities or stochastic processes to operate, requires that iV >> I. 
However, for particles emitting 7 rays of energy 300 MeV, < I, thus requiring special 
acceleration mechanisms—if the emitting electrons are indeed accelerated at the shock. 


One mor e constraint comes from the polarization of the knot. Most recent polarimetry 
performed bv iMoran et al.l (120131) gives a high polarization degree ~ 60%, with position angle 
aligned with the symmetry axis, indicating strongly a toroidal magnetic held . Here we nee d 


to be careful about the depolarization ehect due to relativistic kinematics flLvutikov et ah 
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20031) . What happens then is that Lorentz transformation of the emission of a relativistic 


plasma from the comoving frame to the lab frame prodnces a rotation of the polarization 
vector. For a cnrved emitting snrface, neighboring flnid elements have slightly different 
velocities, thns their polarization vectors experience different amonnts of rotation. As one 
snms the contribntion from the visible snrface (this is essentially what we observe), the 
resnlt is some degree of depolarization. We hnd that for the obliqne shock in the Crab 
Nebnla, nnder a n nltrarelativistic appr oximation, the npper limit of polarization degree is 
similar to that in lLvntikov et al.l (120031) : 56.25% for particle spectral index p = 3 and 43.4% 
for p = 2. Thns, the observed high degree of polarization and its possible variation is 
highly constraining on theoretical models. It wonld be interesting to repeat this observation, 
especially with adaptive optics. 


7. Conclusions 

We have (1) introduced a new approach to analyzing time series of images; (2) discovered 
that key properties of the knot (radial width, tangential width, flnx) are correlated with the 
time-variable separation between the knot and the pnlsar; (3) shown that, with the available 
data, it is not possible to determine a strong correlation between the knot separation and 
the occurrence of y-ray flares; (4) set an upper limit to the low-energy X-ray flux from the 
inner knot; (5) discussed the implications of our observations to set constraints on particular 
elements of the standard shock modeling of the relativistic outflow from the pulsar. These 
include inferences as to wind magnetization a, shock shape parameters such as incident 
angle hi and poloidal radius of curvature Rc, as well as the IR/optical emitting particle 
enthalpy fraction. We found that while the standard shock model gives good agreements 
with observations in many aspects, there remain two puzzles: (a) the angular size of the knot 
relative to the pulsar-knot separation is much smaller than expected; (b) the variable, yet 
high degree of polarization is difficult to reconcile with a highly relativistic outflow. We also 
found that the IR/optical flux of the inner knot is marginally consistent with the scenario 
that the shock accelerates most of the optical emitting particles in the nebula. 
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A. Appendix - Data Analysis 
A.l. Keck 

To study the knot using the high-resolution Keck adaptive-optics (AO) data, we char¬ 
acterize the region around the pulsar using an empirically measured point spread function 
and an analytical model of the knot geometry. 


A. 1.1. PSF Subtraction 

We subtracted the PSF of the pulsar using a nearby PSF star (Figure HOj). During our 
observations, the laser guide star, and so the center of the anisoplanatic patch, was aimed at 
the point half-way between the PSF star and the pulsar. We £t the PSF star, a background 
plus a power-law based seeing-disk to the pulsar. For the £t, we excluded the region in the 
pulsar subimage that surrounds the knot, and used the other 270°. The background was fit 
with a sloping plane model, to allow for the variable nebular background in the pulsar and 
PSF star subimages. The htting was performed with a Levenburg-Marquadt htter. 


A. 1.2. Deconvolution 

We chose not to deconvolve the knot with the empirically measured point spread func¬ 
tion. Tests on several images showed that deconvolution of the shape of the knot made 
insignihcant differences to the fit values and uncertainty. 
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A. 1.3. Knot Model 

To fit the geometric parameters of the knot, we masked out the center diffraction-limited 
core of the pulsar, and used a Levenburg-Marquadt htting technique to £t a two dimensional 
Gaussian to the shape of the Knot. 

We used a model of the form 

^(r,'0) = Sb + Skew{-\{{r - rof / 

We also allowed the center of the coordinate system (xo,|/o) to vary slightly (0.05"). 
To determine the £t uncertainties, we used the 1-a errors from the self-covariance of the fit 
parameters. Results were tabulated in Table 01 


A.2. HST 


A.2.1. HST - Traditional Analysis 


In order to measure the knot properties (position , flux, tangential and radial width) we 
used the SExtractor package flBertin fc Arnontslll990l ). SExtractor has been extensively used 
for the analysis of H ST data, in particula r for the ACS camera (e.g., for the Hubble Ultra 
Deep Field project, iBeckwith et al.l l2006l) . For source detection, we require a minimum 
of hve contiguous pixels with a detection threshold 5a above the root-mean-square (rms) 
background, with a total of 32 deblending subthresholds, and with a contrast parameter of 
0.005, setting the background mesh size to 16x16 ACS pixels. 

The position of a source is evaluated by SExtractor as the barycenter of the source 
brightness distribution. Flu x is c o mput ed within an elliptical aperture, using an imple¬ 
mentation of the method bv iKronI (119801) . Note that the pulsar is not saturated and that 
diffraction spikes from the pulsar that might cross the inner knot only marginally affect its 
flux measurements. The parameters of this ellipse (semiaxes and orientation), evaluated 
using the second moments of the object’s brightness distribution also yield a measure of the 
object morphology. Indeed, for the case of the knot, the direction of the minor axis of the 
ellipse turned out to be consistent with the pulsar-knot direction in all images. Thus, the 
minor axis and the major axis of the ellipse measure the RMS full widths (FWRMS) in the 
radial and tangential directions, respectively. 


In order to assess systematic errors, we performed simulations with the ESO/MIDAS 
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softwar^. We added to the ACS images a “synthetic knot”. To generate such an artihcial 
source, we assumed a two-dimensional Gaussian brightness distribution, with the minor axis 
aligned with the true pulsar-knot direction. The synthetic knot was positioned to the NW 
of the pulsar, opposite but along the true pulsar-knot direction, and at an angular distance 
comparable to that of the true knot. We repeated the exercise by varying the flux, position, 
and morphology of the synthetic knot and we estimated the uncertainties in the parameters 
recovered using SExtractor^ 


Table [TU] lists the best-£t values of the model parameters of the inner knot for each of 
the 17 observations based upon this traditional analysis and modeling. Table [11] shows the 
results of a regression analysis on the listed parameters. 


A.2.2. HST Singular Value Decomposition (SVD) 

As the inner knot is within 1" of the Crab pulsar and varies in position and size, we 
also developed special procedures for characterizing it using the central 121 x 121-pixel 
(6.05" X 6.05") images of the 17 HST/ACS observations. The major steps are (1) to remove 
the pulsar from each image, (2) to use the residual image to generate an image of the inner 
knot, and (3) to characterize the properties of the inner knot. 

Remove pulsar and generate an image. 

In order to remove the pulsar from each of the 17 central 121 x 121-pixel (6.05" x 6.05") 
HST images, we hrst determine the PSF of the HST/ACS using 121 x 121-pixel images 
around each of 19 isolated stars in the held. In doing this, we account for the fact that some 
observations were at a roll angle ~ 180 deg opposite that of the others. After subtracting a 
htted linear gradient from the image of each fFigure fTTl) of the 19 stars and for each of the 
17 observations, we register the resulting 323 = 19 x 17 star images and use singular-value 
decomposition (SVD) to generate a linear basis (lowest term shown in Figure [TT]) describing 
the PSF. The model pulsar image fFigure [TT]) uses the hrst 72 components of the SVD basis. 
After subtracting a linear gradient from each of the 17 central (pulsar) images (Figure [T2|), we 
ht each using the PSF basis to generate 17 residual (pulsar-removed) images (also Figure [T^ . 
each now dominated by the inner knot. 

Characterize inner-knot properties. 


^ https://www.eso.org/sci/software/esomidas/ 
^http://www.astromatic.net/software/sextractor. 
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The measured extrinsic properties of the inner knot are the projected radial separation 
vq and polar angle i/jq of the peak surface brightness with respect to the pulsar. To facilitate 
characterization of the intrinsic properties of the inner knot, we map each inner-knot sub¬ 
image (Figure [12]) onto a ifj-r grid [Figure [T5]) . For an initial model of the surface brightness 
(intensity) distribution of the inner knot, we use a simple bivariate normal distribution: 
S{r, 'll)) = Sk exp(-i((r - ro)Vf^r + (V' - V'o)Vf^p))- 

The major intrinsic properties of the inner knot are its radial dispersion (FWHM^ = 
2.35(Tr for a Gaussian prohle), azimuthal dispersion (or tangential dispersion at = toct.^), 
and peak intensity (surface brightness) Sk- Integrating the intensity over solid angle gives 
the flux Fk = 2naratSk- In htting this model to the 'ip-r image of the inner knot, we include 
a constant-surface-brightness background Sb as a model parameter. 

While this simple model provides adequate estimates to characterize the primary prop¬ 
erties of the knot (Table [6|), as a check we also introduced a somewhat more complicated 
model to deal with minor asymmetries in the 'ij)-r plane. In particular, we took the ra¬ 
dial location and radial width to be weak functions of the azimuthal angle, which we ex¬ 
pand as a Taylor series to second order: r r{'ij)) ro -|- Tq^'i/) — 'ipo) + and 

ar ~ (Jr {'ll)) ~ (Jro + (J^g(V' ~ V'o) + (V' ~ V'o)^- However, the derivatives proved not to alter 

significantly the results. 

Comparison between methods 

The comparison between methods and the reasons for somewhat emphasizing the SVD- 
based results considers the various measured parameters individually. Figure [I4| compares 
measurements of similar variables between the two methods and illustrates the fact that the 
two methods. From the hgure we see that, with the possible exception of the knot flux, both 
methods yield measurements that track each other apart from a scale factor. This leads to 
the result that the pulsar-knot separation determines all the intrinsic properties of the knot. 

Differences for tq are quite small, whereas the knot flux and size estimates differ by 
factors of ~ 2. We would expect differences in the size estimates as the two analyses are 
using different measures of knot size—the hrst, RMS full width (FWRMS); the second, 
FWHM—which are not simply related for a non-Gaussian prohle. That the huxes differ 
is not perhaps so surprising as its measnrement is sensitive to the method of background 
subtraction and the impact of emission in the wings of the inner knot’s prohle. A signihcant 
diherence between the two methods is the smaller uncertainty in tq ahorded by the SVD 
approach, which removes those SVD components that represent noise, thus enhancing signal 
to noise. 
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A. 3. Chandra 

Figure |5] shows the summed image at pulse minimum from the 6 “flaring” observations 
(ObsIDs 14684-16247 in Table [3|). Prior to the binning used to produce the hgure, we extract 
the counts in the 3 regions shown. It is difficult, if not impossible, to make use of the absolute 
number of counts as the different instrument conhgurations are all uncalibrated. However, we 
can make use of the relative numbers of counts from the different regions. Thus, the central 
region serves as an indicator of the number of pulsar counts which can then be compared 
with the other observations to derive physical fluxes (assuming of course that the phase 
averaged pulsar flux has not varied). The region to the southeast provides the upper limit 
for the flux in the knot and we use the region to the northwest to estimate the background, 
comprised mainly of that due to the wings of the pulsar PSF. 

The alert reader will notice a slight excess of counts roughly 0.6" to the southwest of the 
pulsar. This feature is seen both in the phase-averaged pulse and in the data from ObsID 
11245. (The asymmetric effect of the blade inserted for Obs ID 9765 does not allow these 
data to be used in making this comparison.) Therefore, we conclude that this feature is part 
of the pulsar PSF. A more subtle effect is also present as there is a slight asymmetry in the 
PSF that causes the region we associate with the inner knot to be somewhat brighter than 
the background during observations taken in the fall. In the spring, when the spacecraft 
has (naturally) been set at a roll angle that differs by ~ 180 deg from that in the fall, the 
background region is brighter than that we associate with the location of the inner knot. 
This effect is mainly seen in the phase averaged data, and by comparing data from the ObsID 
9765 (spring) and ObsID 11245 (fall) observations. Both effects described in this paragraph 
are relatively minor, are mainly due to the fact we are working slightly below the spatial 
resolution of Chandra, and have thus been ignored in setting the upper limits in Table [S] 
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Table 1: Summary of the Keck observations 


# 

Date 

MJD“ 

Instrument^ 

Exposure (s) 

K1 

2012-02-08 

55965.285 

NIRC2-wide 

420 

K2 

2012-03-05 

55991.319 

NIRC2-wide 

420 

K3 

2012-12-22 

56283.408 

NIRC2-wide 

1250 

K4 

2012-12-23 

56284.410 

NIRC2-narrow 

870 

K5 

2012-12-24 

56285.388 

NIRC2-narrow 

1350 

K6 

2012-12-25 

56286.390 

NIRC2-wide 

1125 

K7 

2013-02-06 

56329.308 

NIRC2-narrow 

2020 

K8 

2013-10-22 

56587.448 

NIRC2-narrow 

640 

K9 

2014-01-09 

56666.274 

NIRC2-wide 

1560 

KIO 

2014-01-09 

56666.377 

NIRC2-narrow 

2580 

Kll 

2014-01-17 

56674.325 

NIRC2-wide 

2270 

K12 

2014-01-17 

56674.394 

NIRC2-narrow 

1260 


“ Quoted MJD is at the midpoint of the observation. 

^ Plate scale is 0.04"/pixel for NIRC2-wide; 0.01"/pixel for NIRC2-narrow. 
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Table 2: Summary of the HST observations 


# 

Date 

MJD“'’ 

HI 

2012-01-08 

55934.797 

H2 

2012-02-10 

55967.170 

H3 

2012-03-12 

55998.414 

H4 

2012-04-22 

56039.850 

H5 

2012-08-16 

56155.336 

H6 

2012-09-10 

56180.960 

H7 

2013-01-10 

56302.992 

H8 

2013-02-24 

56347.055 

H9 

2013-03-06 

56357.023 

HIO 

2013-04-01 

56383.320 

Hll 

2013-04-14 

56396.113 

H12 

2013-08-13 

56517.710 

H13 

2013-10-20 

56585.793 

H14 

2013-10-29 

56594.700 

H15 

2013-12-01 

56627.273 

H16 

2014-01-20 

56677.530 

H17 

2014-04-13 

56760.900 


“ Quoted MJD is at the midpoint of the observation. 
^ Plate scale is 0.05"/pixeh 


Table 3: Summary of the Chandra observations 


# 

ObsID 

Date 

MJD“ 

Exposure (s) 

Cl 

09765 

2008-01-22 

54487.670 

95082 

C2 

11245 

2010-11-16 

55516.294 

22039 

C3 

14684 

2013-03-05 

56356.972 

19485 

C4 

14686 

2013-03-10 

56361.338 

20052 

C5 

14687 

2013-03-17 

56368.324 

16718 

C6 

16244 

2013-10-19 

56584.658 

20048 

C7 

16246 

2013-10-22 

56587.515 

20013 

C8 

16247 

2013-10-24 

56589.007 

19770 


“ Quoted MJD is at the start of the observation. 
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Table 4: Knot properties based upon analysis of the Keck data 



V'o 

ro 

FWHM^ 

FWHMt 

Fk 

# 

o 

// 

ff 

// 

Fp 

K1 

118.50 ±0.16 

0.6711 ±0.0015 

0.3931 ± 0.0040 

0.4719 ±0.0051 

0.0639 ± 0.0045 

K2 

121.00 ±0.13 

0.6577 ±0.0015 

0.4121 ±0.0035 

0.4252 ± 0.0042 

0.0626 ± 0.0067 

K3 

122.64 ±0.11 

0.6893 ±0.0011 

0.3215 ±0.0028 

0.3948 ± 0.0039 

0.0804 ± 0.0074 

K4 

120.96 ±0.18 

0.7118 ±0.0021 

0.3401 ±0.0051 

0.4081 ± 0.0059 

0.0473 ± 0.0003 

K5 

121.20 ±0.07 

0.6981 ±0.0007 

0.2983 ±0.0019 

0.3767 ±0.0024 

0.0555 ± 0.0005 

K6 

120.79 ±0.12 

0.7164 ±0.0010 

0.3375 ± 0.0029 

0.4238 ± 0.0042 

0.0599 ± 0.0085 

K7 

121.93 ±0.12 

0.7756 ± 0.0013 

0.3355 ± 0.0034 

0.4049 ± 0.0038 

0.0662 ± 0.0004 

K8 

120.33 ±0.06 

0.5722 ± 0.0004 

0.2216 ±0.0011 

0.3172 ±0.0015 

0.0621 ±0.0005 

K9 

121.46 ±0.08 

0.6543 ± 0.0006 

0.2958 ±0.0018 

0.4030 ± 0.0024 

0.0708 ± 0.0047 

KIO 

115.77 ±0.03 

0.6257 ±0.0003 

0.2689 ±0.0007 

0.3683 ± 0.0009 

0.0663 ± 0.0006 

Kll 

121.49 ±0.10 

0.6371 ± 0.0008 

0.3119 ±0.0022 

0.3787 ±0.0029 

0.0647 ±0.0078 

K12 

118.99 ±0.09 

0.6438 ± 0.0009 

0.3248 ± 0.0023 

0.3766 ± 0.0025 

0.0574 ±0.0001 
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Table 5: Results of regression analyses for the 12 Keck observations 


Power-law for ri = 0.6692" 
Property 

Unit 

pin) 

Power-law index q 

Probability 

FWHM^ 

n 

0.3180 ±0.0132 

1.19 ±0.54 

5.4E-2 

FWHMi 

n 

0.3941 ± 0.0093 

0.75 ±0.31 

3.8E-2 

Sk 

Fp/n 

0.439 ± 0.033 

-2.14 ±0.97 

5.2E-2 

Fk 

Fp 

0.0626 ± 0.0025 

-0.20 ±0.53 

7.1E-1 

Linear for < tq >= 0.6711" 
Property 


'ipo{< ro >) 

d'lpo/dvQ 


ipo 

7 " 

120.42 ±0.51 

15.2 ± 10.2 

1.7E-1 


Table 6; Knot properties based upon SVD analysis of the 17 HST observations (see Appendix) 


# 

4’o 

O 



^0 

// 


FWHMr 

FWHMt 

// 

Sk 

e/s/(")^ 

Fk 

e/s 


F 

e, 

p 

/s 


Fk 

Fp 


HI 

120.54 


0.40 

0.6819 

± 

0.0023 

0.1526 


0.0056 

0.3221 


0.0116 

1064 ± 33 

59.8 


3.6 

1425 

± 82 

0.0420 

± 

0.0035 

H2 

119.45 


0.41 

0.6535 


0.0024 

0.1536 


0.0058 

0.3071 


0.0114 

1008 ± 32 

53.8 


3.4 

1433 

± 80 

0.0375 


0.0032 

H3 

120.37 


0.33 

0.6708 


0.0022 

0.1751 


0.0052 

0.3121 


0.0092 

902 ± 23 

55.0 


2.6 

1471 

± 64 

0.0374 


0.0024 

H4 

121.74 


0.41 

0.6356 


0.0024 

0.1465 


0.0056 

0.2942 


0.0110 

1053 ± 34 

52.3 


3.2 

1276 

± 80 

0.0410 


0.0035 

H5 

117.90 


0.39 

0.6166 


0.0020 

0.1406 


0.0048 

0.3043 


0.0102 

1256 ± 36 

61.1 


3.4 

1231 

± 84 

0.0497 


0.0044 

H6 

117.62 


0.41 

0.6577 


0.0024 

0.1456 


0.0058 

0.3029 


0.0112 

1102 ± 36 

56.7 


3.5 

1184 

± 86 

0.0479 


0.0046 

H7 

121.17 


0.59 

0.7240 


0.0035 

0.1704 


0.0084 

0.3539 


0.0176 

675 ± 28 

44.7 


3.9 

1422 

± 77 

0.0314 


0.0032 

H8 

119.91 


0.50 

0.7516 


0.0031 

0.1751 


0.0076 

0.3533 


0.0158 

667 ± 25 

44.5 


3.4 

1421 

± 65 

0.0313 


0.0028 

H9 

120.88 

± 

0.52 

0.7491 


0.0033 

0.1740 


0.0080 

0.3362 


0.0166 

728 ± 29 

45.8 


3.8 

1359 

± 73 

0.0337 


0.0033 

HIO 

120.60 

± 

0.32 

0.6614 


0.0019 

0.1598 


0.0046 

0.3139 

± 

0.0090 

1098 ± 27 

63.3 


3.0 

1295 

± 72 

0.0489 


0.0036 

Hll 

119.39 

± 

0.36 

0.6557 


0.0021 

0.1639 


0.0050 

0.3220 


0.0098 

1090 ± 28 

64.9 


3.3 

1509 

± 77 

0.0430 


0.0031 

H12 

119.28 

± 

0.41 

0.6324 


0.0023 

0.1691 


0.0056 

0.3299 


0.0110 

1028 ± 29 

63.4 


3.6 

1283 

± 83 

0.0494 


0.0042 

H13 

117.85 

± 

0.36 

0.5687 

± 

0.0019 

0.1426 


0.0046 

0.2739 


0.0082 

1668 ± 44 

75.6 


3.8 

1222 ± 106 

0.0618 


0.0062 

H14 

118.82 

± 

0.41 

0.5849 

zh 

0.0022 

0.1441 


0.0052 

0.2886 

± 

0.0098 

1583 ± 48 

75.6 


4.4 

1285 ± 118 

0.0589 


0.0064 

H15 

118.59 

± 

0.52 

0.6624 

± 

0.0026 

0.1476 


0.0061 

0.3393 

± 

0.0144 

1105 ± 40 

61.8 


4.4 

1353 ± 100 

0.0457 


0.0047 

H16 

118.82 

± 

0.41 

0.6299 

± 

0.0026 

0.1640 


0.0064 

0.2961 


0.0106 

1084 ± 35 

62.1 


3.7 

1419 ± 98 

0.0438 


0.0040 

H17 

120.88 

± 

0.41 

0.6448 

ih 

0.0024 

0.1611 


0.0058 

0.3212 


0.0110 

1164 ± 35 

69.6 


4.0 

1177 ± 99 

0.0591 


0.0060 
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Table 7: Results of regression analyses for the 17 HST observations based on the data in 
Table [6] 


Power-law for ri = 0.6559" 
Property 

Unit 

p(ri) 

Power-law index q 

Probability 

FWHM^ 

// 

0.1575 ±0.0022 

0.72 ±0.19 

1.7E-3 

FWHMt 

n 

0.3152 ±0.0029 

0.80 ±0.13 

1.4E-5 

Sk 

F.m 

0.782 ±0.023 

-3.66 ±0.40 

1.7E-7 

Fk 

Fp 

0.0440 ± 0.0012 

-2.39 ±0.37 

9.8E-6 

Linear for < tq >= 0.6577" 
Property 


'0o(< >) 

d'ljjo/dro 


^0 

7" 

119.64 ±0.27 

13.5 ±5.5 

2.7E-2 


Table 8: Results of the Chandra analysis 




Pulse Average 

Pulse Minimum 

Fk/Fp 

ObsID 

Pulsar 

Knot 

Background 

Pulsar 

Knot 

Background 

3a 

09765 

51844 

2963 

3704 

845 

57 

72 

< 0.0022 

11245 

5630 

1018 

851 

88 

16 

12 

< 0.0094 

last 6 

13549 

2622 

2503 

203 

39 

44 

< 0.0067 
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Table 9: Fermi/LAT 12-hr-average fluxes and pulsar-knot separations as measured with Keck 
and HST (SVD analysis) 


# 

Date 

^mid 

MJD 

ro 

// 

Fermi/LAT >100 MeV 
10“® ph/s/cnF 

HI 

2012-01-08 

55934.797 

0.6819 ±0.0023 

2.45 

(±0.47, -0.43) 

K1 

2012-02-08 

55965.280 

0.6711 ±0.0015 

2.58 

(±0.60, -0.54) 

H2 

2012-02-10 

55967.171 

0.6535 ± 0.0024 

2.79 

(±0.65, -0.58) 

K2 

2012-03-05 

55991.311 

0.6577 ±0.0015 

2.78 

(±0.46, -0.43) 

H3 

2012-03-12 

55998.414 

0.6708 ± 0.0022 

3.55 

(±0.52, -0.49) 

H4 

2012-04-22 

56039.850 

0.6356 ± 0.0024 

2.18 

(±0.56, -0.50) 

H5 

2012-08-16 

56155.336 

0.6166 ±0.0020 

3.25 

(±0.56, -0.52) 

H6 

2012-09-10 

56180.960 

0.6577 ±0.0024 

3.86 

(±0.50, -0.47) 

K3 

2012-12-22 

56283.381 

0.6893 ±0.0011 

3.13 

(±0.68, -0.63) 

K4 

2012-12-23 

56284.403 

0.7118 ±0.0021 

3.51 

(±0.67, -0.62) 

K5 

2012-12-24 

56285.379 

0.6981 ±0.0007 

3.29 

(±0.67, -0.62) 

K6 

2012-12-25 

56286.370 

0.7164 ±0.0010 

2.74 

(±0.62, -0.56) 

H7 

2013-01-10 

56302.992 

0.7240 ± 0.0035 

3.22 

(±0.66, -0.66) 

K7 

2013-02-06 

56329.248 

0.7756 ±0.0013 

3.93 

(±0.75, -0.70) 

H8 

2013-02-24 

56347.055 

0.7516 ±0.0031 

5.44 

(±0.72, -0.68) 

H9 

2013-03-06 

56357.023 

0.7491 ± 0.0033 

10.09 (±0.52, -0.51) 

HIO 

2013-04-01 

56383.320 

0.6614 ±0.0019 

3.50 

(±0.68, -0.63) 

Hll 

2013-04-14 

56396.113 

0.6557 ±0.0021 

2.89 

(±0.59, -0.54) 

H12 

2013-08-13 

56517.710 

0.6324 ±0.0023 

2.65 

(±0.46, -0.43) 

H13 

2013-10-20 

56585.793 

0.5687 ±0.0019 

3.71 

(±0.26, -0.25) 

K8 

2013-10-22 

56587.443 

0.5722 ± 0.0004 

1.25 

(±1.73, -1.73) 

H14 

2013-10-29 

56594.700 

0.5849 ± 0.0022 

7.46 

(±1.16, -1.07) 

H15 

2013-12-01 

56627.273 

0.6624 ±0.0026 

2.15 

(±0.53, -0.47) 

K9 

2014-01-09 

56666.316 

0.6543 ± 0.0006 

3.38 

(±0.53, -0.50) 

KIO 

2014-01-09 

56666.438 

0.6257 ±0.0003 

3.54 

(±0.51, -0.48) 

Kll 

2014-01-17 

56674.356 

0.6371 ±0.0008 

3.84 

(±0.53, -0.50) 

K12 

2014-01-17 

56674.400 

0.6438 ± 0.0009 

3.47 

(±0.44, -0.42) 

H16 

2014-01-20 

56677.530 

0.6299 ± 0.0026 

2.85 

(±0.43, -0.41) 

H17 

2014-04-13 

56760.900 

0.6448 ± 0.0024 

3.96 

(±0.78, -0.71) 
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Table 10: Results of the Traditional HST Analysis 


# ro FWRMS, FWRMSt Fk 




If 



ff 



ff 


e/s 


HI 

0.660 

± 

0.012 

0.287 

± 

0.019 

0.551 

± 

0.037 

116.5 

± 

2.9 

H2 

0.634 

± 

0.012 

0.283 

± 

0.020 

0.538 

± 

0.038 

110.9 

± 

2.8 

H3 

0.649 

± 

0.013 

0.333 

± 

0.021 

0.594 

± 

0.038 

131.6 

± 

3.3 

H4 

0.615 

± 

0.012 

0.255 

± 

0.020 

0.489 

± 

0.038 

90.4: 

±: 

2.3 

H5 

0.588 

± 

0.012 

0.278 

± 

0.019 

0.560 

± 

0.038 

112.9 

± 

2.8 

H6 

0.635 

± 

0.012 

0.290 

± 

0.020 

0.552 

± 

0.037 

111.6 

± 

2.8 

H7 

0.706 

± 

0.010 

0.288 

± 

0.021 

0.589 

± 

0.044 

97.6: 

i: 

2.4 

H8 

0.734 

± 

0.017 

0.336 

± 

0.024 

0.658 

± 

0.046 

109.2 

± 

2.7 

H9 

0.732 

± 

0.013 

0.319 

± 

0.022 

0.650 

± 

0.046 

108.1 

± 

2.7 

HIO 

0.641 

± 

0.009 

0.299 

± 

0.020 

0.580 

± 

0.039 

117.6 

± 

2.9 

Hll 

0.637 

± 

0.011 

0.314 

± 

0.020 

0.561 

± 

0.036 

129.2 

± 

3.2 

H12 

0.613 

± 

0.012 

0.258 

± 

0.017 

0.578 

± 

0.039 

121.3 

± 

3.0 

H13 

0.547 

± 

0.009 

0.224 

± 

0.015 

0.491 

± 

0.033 

119.2 

± 

3.0 

H14 

0.552 

± 

0.009 

0.236 

± 

0.015 

0.511 

± 

0.033 

124.8 

± 

3.1 

H15 

0.631 

± 

0.012 

0.268 

± 

0.018 

0.588 

± 

0.040 

115.2 

± 

2.9 

H16 

0.614 

± 

0.011 

0.266 

± 

0.025 

0.519 

± 

0.049 

116.9 

± 

2.9 

H17 

0.624 

± 

0.012 

0.259 

± 

0.024 

0.549 

± 

0.050 

127.8 

± 

3.2 


Table 11: Results of regression analyses for the 17 HST observations based on the traditional 
data analysis 


Power-law for ri = 0.6342" 
Property 

Unit 

p(ri) 

Power-law index q 

Probability 

FWRMS, 

ff 

0.2801 ±0.0046 

1.14 ±0.21 

6.0E-5 

FWRMSt 

ff 

0.5605 ±0.0067 

0.85 ±0.15 

4.9E-5 

Fk 

e/s 

114.8 ±2.6 

-0.39 ±0.29 

1.9E-1 
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Fig. 1.— Individual images of the inner knot, with the pulsar removed, from Keck NIRC2 
with laser-guide-star adaptive optics (AO). In each image, the PSF of the pulsar has been 
subtracted using the scaled PSF of a nearby star and replaced with a gray circle at the 
original location of the pulsar. Each image is 2.2" by 2.2". Images are from the NIRC2 wide 
camera (0.04" pixels) or the NIRC2 narrow camera (0.01" pixels). All images were taken 
using the K' hlter. A red dot marks the average position of the inner knot across all 12 
observations. 
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2014 - 04-13 



2012 - 04-22 



Fig. 2.— HST-SVD-Processed images of the Crab pulsar and inner knot after removing 
effects of the HST/ACS WFC point spread function. A red dot marks the average position 
of the inner knot across all 17 observations. 
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Fig. 3.— Variation of the Keck-measured knot properties versus separation of the inner knot 
from the pulsar. 
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Fig. 4.— Variation of the HST-SVD-measured properties with separation of the inner knot 
from the pulsar. 
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A RA (arcsec) 


Fig. 5.— Chandra image at pulse minimum from the sum of the last 6 flaring data sets in 
Table |3l The 3 small black circles show 0.3"-radius extraction regions, the central of which 
is used to estimate the pulsar flux. The circle, at a position angle ijj = 120° east of north 
and 0.65" from the pulsar, is roughly centered on the average optical location of the knot. 
Data in the opposite small circle (the one to the NW where N is up) are used to estimate 
background. Large circles show the best-£t Gaussian to the phase-average pulsar data at 
intensity levels of 10, 20, 40, 80, 160, 320 and 640 cts/pixel illustrating the level of impact of 
the PSF at the site of the knot—i.e., slight but non-negligible. The intensity level is shown 
by the grayscale bar on the right. 
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Fig. 6.— Fermi/LAT 12-hr average flux (upper panel) and pulsar-knot separation (lower 
panel) versus time, for Keck (red) and HST (green) observations. 
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Fig. 7 .— Fermi/LAT 12-hr average fluxes during Keck (red) and HST (green) observations 
at the time of the largest 7 -ray flares in 2013 March (upper panel) and in 2013 October 
(lower panel). 
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Fig. 8.— Fermi/LAT 12-hr-average fluxes centered on the times of Keck (red) and HST 
(green) measurements of pulsar-knot separations. 
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Plane of the slcy 



Fig. 9.— Geometry of the termination shock on the meridional plane determined by the 
pulsar rotation axis and the observer’s line of sight toward the pulsar. The pulsar is located 
at point P. The angle b etween line of s i ght an d pulsar rotation axis is 6*ob = 27r/3, consistent 
with that measured bv iNg fc Romani (120081) . The knot is projected on the plane of the sky 
at an angle tq ^ 0.65" southeast of the pulsar (represented by point K here). The angle 
between the upstream velocity and the shock surface is (5i. The outflow is deflected from the 
radial direction by an angle A. 
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Pulsar & Knot PRF Star 



K8: 2013-10-22 NIRC2-narrow 


Fig. 10.— The Pulsar (left) and a nearby star (right) from Keck data taken on 2013 October 
22. (a) An H hlter image of the pulsar and knot (in the lower-left) as seen by the N1RC2- 
narrow camera (0.01" pixels and a 10" x 10" held of view). The structure of the knot is 
resolved and separated from the pulsar by a statistically signihcant valley, (b) A nearby (5") 
comparison star, used to establish the point-response-function. Both images are shown with 
a logarithmic color stretch. 
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Fig. 11.— Images illustrating the steps in determining the HST/ACS PSF for this analysis: 
(a) Extract 323 (19 x 17) 121 x 121-pixel (6.05" x 6.05") images and subtract a linear-gradient 
background for 19 isolated stars in these 17 HST/ACS observations of the Crab; (b) register 
images of stars and use singular value decomposition (SVD) to determine a basis describing 
the PSF (Erst SVD component shown); and (c) model the pulsar image using the hrst 72 
SVD components. 
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Fig. 12.— Images further illustrating steps in processing HST/ACS images of the Crab pulsar 
and inner knot: (a) Subtract a linear-gradient background in central 121 x 121 pixel (6.05" x 
6.05") image of pulsar and inner knot; (b) remove pulsar using SVD model of PSF leaving 
inner knot; (c) remove inner knot using its SVD model leaving residual background; (d) 
add SVD model of inner knot to residual background; (e) apply Richardson-Lucy algorithm 
to generate deconvolved image of inner knot (and residual background); and (f) synthesize 
25 X 25 pixel (1.25" x 1.25") sub-image of registered (unresolved) pulsar and SVD model of 
inner knot. 
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Fig. 13.— Mapping of processed images of the inner knot onto a ip-r grid. Here ifj is the 
angular displacement in radian from the knot’s centroid; r, the radial distance from the 
pulsar in (0.05") pixels. 
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Fig. 14.— Comparison of HST-SVD and HST-traditional measured properties of the inner 
knot. 

































































